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Abstract 

The spinless relativistic Coulomb problem is the bound-state problem for the spinless 
Salpeter equation (a standard approximation to the Bethe-Salpeter formalism as well 
as the most simple generalization of the nonrelativistic Schrodinger formalism towards 
incorporation of relativistic effects) with the Coulomb interaction potential (the static 
limit of the exchange of some massless bosons, as present in unbroken gauge theories). 
The nonlocal nature of the Hamiltonian encountered here, however, renders extremely 
difficult to obtain rigorous analytic statements on the corresponding solutions. In view 
of this rather unsatisfactory state of affairs, we derive (sets of) analytic upper bounds 
on the involved energy eigenvalues. 
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1 Introduction: The Spinless Salpeter Equation 

Maybe one of the most straightforward generalizations of the standard 
nonrelativistic quantum theory towards the reconciliation with all the 
requirements imposed by special relativity is represented by describing 
the quantum systems under consideration by the well-known "spinless 
Salpeter equation." Consider a quantum system the dynamics of which 
is governed by a by assumption self- adjoint Hamiltonian H of the form 

H = T + V : (1) 

where T denotes the square-root operator of the relativistic expression 
for the free (kinetic) energy of some particle of mass m and momentum 
P- 

T — \Jp 2 + m 2 , 

and V = V(x) represents some arbitrary coordinate-dependent static 
interaction potential. The eigenvalue equation for this Hamiltonian H, 

H\xk) = E k \ Xk ) , k = 0,1,2,... , 

for Hilbert-space eigenvectors \xk) corresponding to energy eigenvalues 

E = (Xk\H\xk) 
{Xk\Xk) 

is nothing else but the one-particle spinless Salpeter equation. Because 
of the nonlocality of this operator H, that is, more precisely, of either 
the kinetic-energy operator T in configuration space or the interaction- 
potential operator V in momentum space, it is hard to obtain analytic 
statements from this equation of motion. (For the "translation" of the 
equal-mass two-particle problem to the one-particle problem discussed 
at present, see Appendix [A].) 

2 The Spinless Relativistic Coulomb Problem 

Of particular importance in all physics are the (spherically symmetric) 
central potentials, which depend only on the radial coordinate r = |x|. 
Among these, the most prominent one is the Coulomb potential Vc{r) } 
which is parametrized by some (dimensionless) coupling constant a: 



V(x) = Vc(r) = , a > 



(2) 
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The (semi-)relativistic Hamiltonian ([[]) with the Coulomb interaction 
potential Vq in (|2j) defines the "spinless relativistic Coulomb problem." 
In the past, the interest in this spinless relativistic Coulomb problem 
has undergone an eventful history. (For a rather comprehensive review, 
consult Ref. |Tj.) Let us merely sketch in the following some highlights. 

First of all, from an examination [j2j of the spectral properties of the 
operator ([I]), (|2j) one may infer the existence of its Friedrichs extension 
up to the critical value 



a c = - 

7T 



of the involved coupling constant a and read off a lower bound on the 
corresponding ground-state energy eigenvalue Eq, namely, 



Eq > m 

which, later on, has been slightly improved to H 
Eq > m 



1 + VI -4a 2 r 1 

\ ^ for a < - 

\ 2 2 



The analytic solutions for the wave functions of those eigenstates |x) of 
the Hamiltonian ([I]), (|2j) which correspond to vanishing orbital angular 
momentum have been constructed [|j . The attempt in [|J to determine 
simultaneously the respective (set of) exact energy eigenvalues without 
actually solving this spinless Salpeter equation, however, failed [§, |6j]. 
Therefore, as far as analytic statements about the relativistic Coulomb 
problem, in particular, its energy eigenvalues, are concerned, up to now 
one has to content oneself with a few series expansions of these energy 
eigenvalues Ek in powers of the coupling constant a |7], |8j] , which then 
are, of course, only significant for a region of rather small values of a. 



3 Analytic Upper Bounds on Energy Eigenvalues 

Without a closed form of all the energy eigenvalues Ek of the spinless 
relativistic Coulomb problem at hand, it is highly desirable to have, at 
least, analytic expressions for upper bounds on these at one's disposal, 
in order to estimate the reliability of approximative solutions or series 
expansions. 

The theoretical basis as well as the primary tool for the derivation of 
rigorous upper bounds on the eigenvalues of some self- adjoint operator 



is, beyond doubt, the so-called "min-max principle" Q9j]. An immediate 
consequence of this min-max principle is the following statement: Let 
H be a self-adjoint operator that is bounded from below [as, according 
to Eq. (|3|), evidently holds for the (semi-)relativistic Hamiltonian in (JTJ) 
with a Coulomb-type interaction potential (@)]. Let Ek, k = 0, 1, 2, . . ., 
denote the eigenvalues of H, ordered according to Eq < E\ < Ei < . . .. 
Let Dd be some d-dimensional subspace of the domain of H . Then the 
kth eigenvalue (counting multiplicity) of H , Ek, satisfies the inequality 



A peculiarity of the (spinless relativistic) Coulomb problem is that 
there is only one dimensional parameter, namely, the particle mass m. 
As a consequence of this, for a vanishing particle mass, i. e., for m = 0, 
the totality of eigenvalues Ek of the Hamiltonian ([I]) with interaction 
potential (0) collapses, already on dimensional grounds, necessarily to 
Ek = for all k = 0, 1, 2, . . .. This fact is also clearly demonstrated by 
application of the "relativistic virial theorem" proven in Refs. [JT^, |TTf| . 
Accordingly, in the course of an investigation of the spinless relativistic 
Coulomb problem, it is sufficient to focus one's interests to the special 
case of a nonvanishing particle mass m: 



In view of the above, when searching for bounds, our intention must 
be to avoid in some way or other the problematic square-root operator 
in order to deal with more manageable Hamiltonians. 

3.1 The "Schrodinger" bound 

From the positivity of the square of the obviously self-adjoint operator 



E k < sup 




k = 0,1,2,... . 



(4) 



m > . 



T — m 



< (T-m) 2 

= T 2 + m 2 -2mT 
= p 2 + 2m 2 -2mT 



one obtains, for the free (or kinetic) energy T, the operator inequality 
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and thus, for the generic Hamiltonian H in ([I]), the operator inequality 

H < H s , (5) 
where Hs denotes the "Schrodinger" Hamiltonian 

H s = m + ^- + V . (6) 
2 m 

Applying to the energy eigenvalues of the Hamiltonian H in Eq. ([I]) 
first the min-max principle in the form given by Eq. and after that 
the operator inequality (|5j), we find 

{Xk\H\Xk} 
{Xk\Xk) 

\H 



Eh — 



< sup 



< sup 



Now, let us assume that the (k + l)-dimensional subspace D^+i in this 
inequality is spanned by the first k + 1 eigenvectors of the Schrodinger 
Hamiltonian Hs-, that is, by precisely those eigenvectors of H$ which 
correspond to the first k + 1 energy eigenvalues £3,0 > ^s,i? • • • ? ^s,fe when 
all eigenvalues of Hs are ordered according to Es,o < £s,i < Es,2 < ■ ■ ■■ 
In this case, the right-hand side of the above inequality is nothing else 
but the kth "Schrodinger" energy eigenvalue Es 



sup = ^S,fc • 

Consequently, any single energy eigenvalue E\. of the spinless Salpeter 
equation is bounded from above by its "Schrodinger" counterpart i?s,fc : 

Ek < E S} k ■ (7) 

For the Coulomb potential (|2|), the energy eigenvalues required here 
are well known: 

Es - = m t 1 - £) • (8) 

where the total quantum number n is given in terms of both radial and 
orbital angular momentum quantum numbers n T and £, respectively, by 

n = n r + £ + 1 , n r = 0, 1, 2, . . . , I = 0, 1, 2, ... . 
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3.2 A "squared" bound 

One might be tempted to try to find an improvement of the bound ([7]) 
by considering the square of the Hamiltonian H, that is, the operator 

Q = H 2 = T 2 + V 2 + TV + VT . 

The eigenvalue equation for this squared Hamiltonian Q will, of course, 
be solved by the same set of eigenvectors \xk) as the one for the original 
Hamiltonian H with, however, the squares of the corresponding energy 
eigenvalues as the eigenvalues of Q: 

Q\Xk) = E 2 k \ Xk ) , k = 0,1,2,... . 

From the positivity of the square of the self-adjoint operator T — m — V, 

< (T-m-V) 2 

= T 2 + m 2 + V 2 - 2mT + 2mV - TV - VT , 

we obtain, with the help of the obvious relation 

< m < T , 

for the anticommutator TV + V T of kinetic energy T and interaction 
potential V, which appears in the square Q of our Hamiltonian H, the 
operator inequality 

TV + VT < T 2 + m 2 + V 2 - 2mT + 2mV 

- p 2 + 2m 2 + V 2 -2mT + 2mV 

< p 2 + V 2 + 2mV , 

which, when applied to the squared Hamiltonian Q, yields the operator 
inequality 

Q<R, 

with the operator R given by 

R = 2p 2 + m 2 + 2V 2 + 2mV . 

Following the line of argument as given in the preceding subsection, we 
might therefore conclude that the squares of the energy eigenvalues E k 
of the spinless relativistic Coulomb problem are bounded from above 
by the eigenvalues E^k of the operator R: 

El < S Rjk , 
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which entails 




(9) 



Now, only for the Coulomb potential (|2j), this latter operator R is of 
precisely the same structure as the Schrodinger Hamiltonian (0), with, 
however, some kind of "effective" orbital angular momentum quantum 
number L which, for states of (genuine) orbital angular momentum £, 
has to be determined from the relation 



In this case, the eigenvalues £r of the operator R may be easily written 
down: 



with the "effective" total quantum number 

N = n v + L + l , n r = 0,1,2,... . 

Unfortunately, it is easy to convince oneself that the "square" bounds 
(0) obtained in this way lie above and are thus worse than the previous 
Schrodinger bounds. 

3.3 A variational bound 

For all practical purposes, the most efficient manner for the min-max 
principle to come into play is in form of the "Rayleigh-Ritz variational 
technique." Let us illustrate this fact just for the ground-state energy 
eigenvalue Eq. For k = 0, the relation simplifies to 



That is, the ground-state energy Eq is clearly less than or equal to any 
expectation value of the considered Hamiltonian, H . The above upper 
bound may, of course, be optimized by determining the smallest of all 
these expectation values, at least in some chosen Hilbert-space sector. 
Consequently, there is a simple recipe for the derivation of (sometimes 
excellent!) exact upper bounds on this ground-state energy eigenvalue: 



L(L + 1) =£{£+!) + a 2 



the solution of which reads 






minimize the encountered expectation values of the Hamiltonian under 
consideration, 



with respect to a suitably chosen set of Hilbert-space trial vectors |^a), 
which are distinguished from each other by some variational parameter 
A. 

We apply this prescription to our (semi-)relativistic Hamiltonian H . 
However, in order to make life easy, we immediately take advantage of 
a simple inequality for the expectation values of a self-adjoint operator, 
like our kinetic energy T, with respect to arbitrary Hilbert-space states 
in the domain of this operator: 



< 



{4>\T 2 \i/j) 



Employing this inequality, we are able to circumvent the (troublesome) 
expectation values of the square-root operator of the kinetic energy T 
(which, of course, may be always evaluated numerically; see Ref. [|Ij): 



En < 



(^a|^a) 
{j;x\T + V\iJ>x) 



< 



<^a|T 2 |^) (^aI^a) 



{i>\\4>) 



(^a|p 2 |^a 



\ <^a) 



+ m l + 



^x\V\^x) 



For the Coulomb potential (@), the most reasonable choice of trial 
vectors is obviously one for which the coordinate-space representation 
■0a( x ) °f the trial vectors for vanishing radial and orbital angular 
momentum quantum numbers is given by hydrogen-like trial functions: 



^a(x) = exp(— Xr) , A > . 



For these trial functions, the computation of the expectation values in 
the above inequality yields |5f 



s 



Determining the minimum of this latter set of upper bounds, we arrive 
at 

Eq < £var,0 i 

with the variational upper bound for the ground-state energy level Eq 
of the spinless relativistic Coulomb problem given by 

£ V ar,o = rn Vl- a 2 . (10) 

This variational bound, i? vari o, is lower and thus better than the former 
Schrodinger bound on the ground-state energy level, 

E s,o = m (l - yj , 

which may be obtained from Eq. (|8j) for n r = t = 0, implying thereby 
n = 1: 

E v&rfi < E Sfi for a ^ . 

Consequently, the variational technique entails indeed improved upper 
bounds on the energy levels as compared to the Schrodinger estimates. 

3.4 A straightforward generalization 

Our variational upper bound (|TD|) for the ground-state energy level Eq 
can be very easily re-derived and simultaneously extended to arbitrary 
levels of excitation by a generalization of the considerations presented 
in Subsection Introducing an arbitrary real parameter fi (with the 
dimension of mass), we make use of the positivity of the square of the 
obviously self-adjoint operator T — /i, 

< (T-/i) 2 

= T 2 + n 2 -2jiT 

= p 2 + m 2 + /i 2 -2/iT , 

in order to find, for the kinetic energy T, a set of operator inequalities, 

T< p2 + m2 + ^ forall M >0, 
2 fi 

and, consequently, for the (semi-)relativistic Hamiltonian H in Eq. ([T]), 
the set of operator inequalities 



H < Hs(/i) for all /i > , 
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with the Schrodinger- type Hamiltonian H^(fi) given by 

Now, mimicking the line of argument given in Subsection f?7T] involving 
the min-max principle, we conclude that the set of energy eigenvalues 
Ek, k = 0, 1, 2, . . ., of the Hamiltonian iif , if again ordered according to 
Eq < E\ < E2 < . . ., is bounded from above by the (corresponding) set 
of energy eigenvalues £s,fc(/i) of this Schrodinger Hamiltonian Hs(fi), if 
again similarly ordered according to £3,0 (/-0 — Es,i(n) < £3,2 (/-0 < • • 

E k < E s ,k(fj) for all /i > , 

and, consequently, also by the minimum of all these energy eigenvalues: 

E k < mmE s ,M ■ 

For the Coulomb potential (0), these energy eigenvalues Es jn (/i) are 
given by 

m 2 + /i 2 (l 



n 2 



2fi 

with precisely the same total quantum number n as before. Minimizing 
this latter expression with respect to the parameter /i, we finally obtain 



a 2 



mm Es, n (fj) =m\l ^ . 

m>o \ n 2 

For any value of the total quantum number n, these bounds definitely 
improve the Schrodinger bounds (|8|): 

mm E s , n (fi) < Es, n for a ^ . 

Clearly, for /i = m, we immediately recover the Schrodinger approach 
of Subsection p7T| . 



4 Summary 

Horrified by the insight that the discrete spectrum of the Hamiltonian 
consisting of just the relativistic kinetic energy and the static Coulomb 
interaction potential is still not known exactly, we derived by different 
but elementary methods at least a few complete sets of analytic upper 
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bounds on the respective energy eigenvalues. In every individual case, 
the basic idea was to derive an operator inequality for this Hamiltonian 
which guarantees that the expectation values of this Hamiltonian are 
bounded from above by the expectation values of some other operator 
which no longer involves the (problematic) square-root operator of the 
relativistic kinetic energy, and to "embed" this operator inequality into 
the well-known min-max principle for the eigenvalues of a self-adjoint 
operator bounded from below. It should be no great surprise that just 
that method which employs some variational procedure yields the best 
of these bounds. 
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A Equivalence of One-Particle and (Equal-Mass) 
Two-Particle Scenarios 

The two relativistic Coulombic Hamiltonians for one- and two-particle 
problem, 

= VP 2 + m 2 , r = Ixl , 

r 

and 

= 2 VP 2 + M 2 - ^ , i?=|X|, 

R 

respectively, may be equated with the help of a scale transformation as 
follows. Relate the employed phase-space variables, (x, p) and (X, P), 
respectively, by some (in general arbitrary) scale factor A according to 

p = AP, 

X 

X = A ' 

which preserves their fundamental commutation relations: 

[x, p] = [X, P] . 

Fix this scale factor, A, to the particular value A = 2 and identify both 
the mass and the Coulomb coupling strength parameters according to 

m = 2M , 
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You will end up with 
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